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HITTING TIMES FOR GAUSSIAN PROCESSES 

By Laurent Decreusefond and David Nualart 1 

GET /Telecom Paris and University of Kansas 

We establish a general formula for the Laplace transform of the 
hitting times of a Gaussian process. Some consequences are derived, 
and particular cases like the fractional Brownian motion are dis- 
cussed. 

1. Introduction. Consider a zero mean continuous Gaussian process (Xt, 
t > 0), and for any a > 0, we denote by r a the hitting time of the level a 
defined by 

(1.1) T a = inf{t > : X t = a} = M{t >0:X t >a}. 

Thus, the map (o i— > r a ) is left-continuous and increasing, hence, with right 
limits. The map (a i— > r a +) is right continuous where 

r a + = lim T a = inf{t > : Xt > a}. 

bla,b>a 

Little is known about the distribution of r a . It is explicitly known in partic- 
ular cases like the Brownian motion. If X is a fractional Brownian motion 
with Hurst parameter H, there is a result by Molchan [5] which stands that 

P( Ta >t)=t-( 1 - H ^°W 

as t goes to infinity. 

When X is a standard Brownian motion, it is well known that 

(1.2) E(exp(—aT a )) = exp(—aV2a) 

for all a > 0. This result is easily proved using the exponential martingale 

M t = exp(\B t - \\ 2 t). 
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By Doob's optional stopping theorem applied at time t A r a and letting t — ► 
oo, one gets 1 = E(M Ta ) = E(exp(\B Ta — A 2 r a /2)). Since i? Ta = a, we thus 
obtain (1.2). If we consider a general Gaussian process X t , the exponential 
process 

M t = exp(\X t -±\ 2 V t ), 

where V t = E(X?) is no longer a martingale. However, it is equal to 1 plus 
a divergence integral in the sense of Malliavin calculus. The aim of this 
paper is to take advantage of this fact in order to derive a formula for 
.E(exp(— \\ 2 V Ta )). We derive an equation involving this expectation in The- 
orem 3.4, under rather general conditions on the covariance of the process. 
As a consequence, we show that if the partial derivative of the covariance is 
nonnegative, then E(exp(—^X 2 V Ta )) < 1, which implies that V Ta has infinite 
moments of order p for all p > ^ and finite negative moments of all orders. In 
particular, for the fractional Brownian motion with Hurst parameter H > i, 
we have the inequality 

-E(exp —oct 211 ) < exp(— aV2a) 

for all a, a > 0. 

The paper is organized as follows. In Section 2 we present some prelimi- 
naries on Malliavin calculus, and the main results are proved in Section 3. 

2. Preliminaries on Malliavin calculus. Let (Xt,t > 0) be a zero mean 
Gaussian process such that Xq = and with covariance function 

R(s,t)=E(X t X s ). 

We denote by £ the set of step functions on [0, +oo). Let TL be the Hilbert 
space defined as the closure of £■ with respect to the scalar product 

(ilo.tjjip.s])^ = R(t,s). 

The mapping lr 0)t i — > Xt can be extended to an isometry between ri and 
the Gaussian space H\(X) associated with X. We will denote this isometry 
by <p — 

Let S be the set of smooth and cylindrical random variables of the form 
(2.1) F = f(X(<f> 1 ),...,X(<f> n )), 

where n > 1, / S C£°(lR n ) (/ and all its partial derivatives are bounded), and 
<f>i£H. 

The derivative operator D of a smooth and cylindrical random variable 
F of the form (2.1) is defined as the W-valued random variable 

n 8f 

i=l * 



HITTING TIMES FOR GAUSSIAN PROCESSES 3 

The derivative operator D is then a closable operator from L 2 (Q) into 
L 2 (Q;TC). The Sobolev space B 1 ' 2 is the closure of S with respect to the 
norm 

\\F\\l 2 = E(F 2 )+E(\\DF\\ 2 n ). 

The divergence operator 5 is the adjoint of the derivative operator. We 
say that a random variable u in L 2 (Q;7i) belongs to the domain of the 
divergence operator, denoted by Dom 5, if 

\E((DF,u) H )\<c u \\F\\ L 2 {n) 

for any F € S. In this case 5(u) is defined by the duality relationship 

(2.2) E(F5(u)) = E((DF,u) n ), 

for any F £ ED 1 ' 2 . 

Set Vt = R(t, t). For any A > 0, we define 

M t = exp(\X t -±\ 2 V t ). 

Formally, the Ito formula for the divergence integral, proved, for instance, 
in [1], implies that 

(2.3) M t = l + X5(M1 M ), 

where Ml^jj represents the process (s t— ► M s l^ ^(s), s > 0). However, the 
process Ml^ ( i does not belong, in general, to the domain of the divergence 
operator. This happens, for instance, in the following basic example. 

Example 1. Fractional Brownian motion with Hurst parameter H £ 
(0, 1) is a zero mean Gaussian process (B^ 1 ,t > 0) with the covariance 

(2.4) RH{t ,s) = ±(t 2H + s 2H -\t-s\ 2H ). 

In this case, the processes (B^l^ ^ (s),s> 0) and (exp(Ai?^ — ^A 2 s 2 ^)l[ 0) t](s), 
s > 0) do not belong to L 2 (n-,H) if H < \ (see [2]). 

In order to define the divergence of Ml^ t ] and to establish formula (2.3), 
we introduce the following additional property on the covariance function of 
the process X. 

(HO) The covariance function R(t, s) is continuous, the partial derivative 
^j(s, t) exists in the region {0 < s, t, s / t}, and for all T > 0, 



sup 

te[o,r]^o 



8R. . 



ds < oo. 



Notice that this property is satisfied by the covariance (2.4) for all H £ 
(0,1). 
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Define 

(2.5) 6 t M=j(M t -l). 

The following proposition asserts that 5tM satisfies an integration by parts 
formula, and in this sense, it coincides with an extension of the divergence 
of Ml m . 

Proposition 2.1. Suppose that (HO) holds. Then, for any t>0, and for 
any smooth and cylindrical random variable of the form F = f(X tl , . . . , X tn ), 
we have 

( n Bf ( l BR \ 

(2.6) E(F5 t M) = E\T^(X tl ,...,X tn ) ^ M s — (s, U) dsj . 

Proof. First notice that condition (HO) implies that the right-hand side 
of equation (2.6) is well defined. Then, it suffices to show equation (2.6) for 
a function of the form 

f(x u . . . , x n ) = exp X i x ^j ' 

where Aj S M. In this case we have, for all <! t\ <C ■ • ■ < t n , 
^E(F(M t - 1)) 

= ^ ex p|^£^ A i^(^^i)| (exp j^AA;i?(Mi) j - IJ 

11 ft ( I n n \ QR 

= Y1 / e^pl-^XiXjRiti,^) + X^2XiR(s,ti)\Xi — {s,ti)ds 

i=l " I i = l t=l J 

-jC*(£|£<*. *.)«.§<..«)*. 

which completes the proof of the proposition. □ 

In many cases like in Example 1 with H > ^, the process Ml[ 0j t] belongs 
to the space L 2 (f2;7Y), and then the right-hand side of equation (2.6) equals 

E{DF,Ml m ) H . 

In this situation, taking into account the duality formula (2.2), equation (2.6) 
says that Ml[ 0t ] belongs to the domain of the divergence and <5(Ml[ 0t ]) = 
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3. Hitting times. In this section we will assume the following conditions: 

(HI) The partial derivative |jf (s, i) exists and it is continuous in [0, +oo) 2 . 

(H2) limsupj^oo Xt = +oo almost surely. 

(H3) For any < s < t, we have E(\X t - X s \ 2 ) > 0. 

Under these conditions, the process X has a continuous version because 

E(\X t - X s \ 2 ) = R(t,t) + R(s,s) - 2R(s,t) 




< 2\t — s\ sup 

s<u<t 

For any a > 0, we define the hitting time r a by (1.1). We know that 
P{r a < oo) = 1 by condition (H2). Set 

(3.1) S t = sup X s . 

s6[0,t] 

From the results of [6], it follows that, for all t > 0, the random variable St 
belongs to the space D ,2 . Furthermore, condition (H3) allows us to compute 
the derivative of this random variable. 



OR 

du 



(u,t) 



Lemma 3.1. For all t > 0, with probability one, the maximum of the 
process X in the interval [0, t] is attained in a unique point, that is, Ts t = t s + 
and DS t = l[o,r s J- 

Proof. The fact that the maximum is attained in a unique point follows 
from condition (H3) and Lemma 2.6 in Kim and Pollard [4]. The formula 
for the derivative of St follows easily by an approximation argument. □ 

We need the following regularization of the stopping time r a . Suppose 
that ip is a nonnegative smooth function with compact support in (0, +oo) 
and define for any T > 

(3.2) Y= / (f(a)(T a AT) da. 

Jo 

The next result states the differentiability of the random variable Y in the 
sense of Malliavin calculus and provides an explicit formula for its derivative. 

Lemma 3.2. The random variable Y defined in (3.2) belongs to the 
space D > 2 , and 

(3.3) D r Y = -J%(y)l^ Ty] (r)dr y . 
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Proof. Clearly, Y is bounded. On the other hand, for any r > 0, we 
have 

{r a >r} = {S r < a}. 
Therefore, we can write using Fubini's theorem 

Y = j (p(a)yj d6jda = J (j tp(a)dajd6, 

which implies that Y € B 1 ' 2 because Sg £ D 1 ' 2 , and 

D r Y = - [ T \p(Sg)D r Sgd6 = - [ p(S e )l [0>TSo] (r)de. 
Jo Jo 9 

Finally, making the change of variable Sg = y yields 

/•St 

D r Y = - (p(y)l [0>Ty] (r)dT y . n 

J 

Notice that My = exp(AAV — \\ 2 Vy)- Hence, letting t = Y in equation 
(2.5) and taking the mathematical expectation of both members of the equal- 
ity yields 

(3.4) E(M Y ) = l + \E{5 t M\ t=Y ). 

We are going to show the following result which provides a formula for the 
left-hand side of equation (3.4). 

Lemma 3.3. Assume conditions (HI), (H2) and (H3). Then, we have 

( f s T dR \ 

(3.5) E(M Y ) = 1-XE [My J v{v)-q^(X, t v ) dr y \ . 

Proof. The proof will be done in two steps. 

Step 1. We claim that for any function p(x) in Cg°(R) we have 

/ ft rS T QE \ 

(3.6) E(StMp(X)) = -EU M e i/(X)J v(y)^(s,T y )dT y ds). 

We can write Y = J Q ip(Sg) dO, where ip(x) = ip(a) da. Consider an in- 
creasing sequence D n of finite subsets of [0, T] such that their union is dense 
in [0, T] . Set Y n = J T ^(S£) d9, and Sfi = max{X t ,t £ D n n [0, 6]}. Then, Y n 
is a Lipschitz function of {Xt,t £ D n }. Hence, formula (2.6), which holds for 
Lipschitz functions, implies that 

E(5 t Mp(Y n )) = -E(p\Y n ) jf <p(S%) (J* M s ^(s, r s? )ds^ d9^j . 
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The function r — > /q M s ^(s,r) ds is continuous and bounded by condition 
(HI). As a consequence, we can take the limit of the above expression as n 
tends to infinity and we get 

E(8 t Mp(Y)) = -E(p'(Y) J <p(S e )(J* M s ^(s,r Sg ) ds^j d8 

Finally, making the change of variable Sg = y yields (3.6). 
Step 2. We write 



E(S t M\ t=Y ) = e( lim [°° 5 t Mp £ (Y - t) 



dt 



where p £ (x) is an approximation of the identity, and by convention, we 
assume that SfM = if t is negative. We can commute the expectation with 
the above limit by the dominated convergence theorem because 

/oo roo 1 

\8 t M\p £ {Y -t)dt= / -\M t - l\p £ (Y - t) dt 
-oo J —oo A 

<- sup (|M t | + l), 
^ 0<*<T+1 

if the support of p £ (x) is included in [— e, e], and e < 1. Hence, 

/oo 
E(S t Mp £ (Y-t))dt. 
-oo 

Using formula (3.6) yields 

sOStMjfeCr-t)) 

(3.8) ft f f r s T BR \\ 

= -jf E[p' £ {Y-t)M s ^ <p(y) — (s,T y )dT y ))ds. 

Hence, substituting (3.8) into (3.7) and integrating by parts, we obtain 
E(8 t M\ t=Y ) 

= -^jy^ Y - t iJ* M °{C^§ (s,r y )dr y )ds)dt 



ft 00 ( r ST dR 

^ E {j_ 0o P ^ Y ~ t \ Mt J V(y)-^{tlTy)dTy 



dt). 



Notice that 



f ST OR 

Jo ^~Qi:( s > T y^ dT i 



< T sup 

0<s,«<T 



dR 



m 
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Hence, applying the dominated convergence theorem, we get 
E(M Y ) = l + \E(5 t M\ t=Y ) 

= 1 - A lim E^J°° p £ {Y - t) (nf t jJ T <p{y)^{t, T y ) dT y ^j dt 

( P s t dR \ 

= 1 - XE [M Y J <p(y) — (Y, Ty)dT y ). □ 

The next step will be to replace the function ip(x) by an approximation 
of the identity and let T tend to infinity. Notice that (3.5) still holds for 
ip(x) = l[o i 6](x) for any b > 0. In this way we can establish the following 
result. 

Theorem 3.4. Assume conditions (HI), (H2) and (H3). For any a > 
and AgM, we have 

E{M Ty )dy 

\Jo Jo MzT y+ +i * l - z ) T y~ds~( ZT y + + ^ ~ z ^ T y> T ^ dz dT y ) ■ 



(3.9) ( f a f l_ dR 

= a — XE 



Notice that we are not able to differentiate with respect to o, the integral 
in the rightmost expectation of (3.9), because the (random) measure dr y , in 
general, is not absolutely continuous with respect to the Lebesgue measure. 

Proof of Theorem 3.4. Fix a > 0. We first replace the function ip(x) 
by an approximation of the identity of the form <p e (x) = e~ 1 l[o ) i](x/e) in 
formula (3.5). We will make use of the following notation: 

POO 

Y £ ,a= / <p e {x - a){T x AT)dx. 



At the same time we fix a nonnegative smooth function ipix) with compact 
support such that f R ip(a) da = c and we set 



/ E{M Yeta )i>{a)da 



c-X E[ My, 



[ St dR \ 

f tp e (y-a)—(Y £ia ,Ty)dTy\ip(a)da. 



The increasing property of the function x — > r x implies that r a + AT < Y £A < 
T a +e A T. Hence, Y £ converges to r a + AT as e tends to zero. Thus, almost 
surely, we have 

\imM Ye = exp(AX r +AT - \X 2 V T +AT ). 



£^0 
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\xm[_E(M Ye , a )il>(a)da= f £(exp(AX r AT - ±\ 2 V Ta+AT )Ma) da. 



Now, set F(t) = M t ^(t,T y ). Then, assuming that <p e (x) = e~ 1 l[ ^(x/e), 
we have 

[V OR 

J Ve(y-a)My eia -^(Y eia ,T v )il>(a)da 

ra+e 



-J" i 

£ J y—s 



F 



[0,1] 



y-g 

£ 



F 



-[0,1] 



x — a 



A T) dx ) Y>(a) da 



Ty+e^-eTj A T) d£j lj}(y - £7]) dl] 

jf F^J\r y+e ^ ev AT)dC + J (r y+e ^ ev AT)d^(y-er ] )dr ] . 



'V 

As e tends to zero, this expression clearly converges to 

■i 



*l>(y) / F(r,(r y A T) + (1 - V )(r y+ A T)) dr,. 
Jo 

So, we have proved that 



OR 

Km / My ^(y - o)— (y 6|(1I t,)^(o) da 



(3.10) 



WVy) J M zry++{1 _ z)Ty —{ZT y ++{l-Z)Ty,Ty)dZ. 



In order to complete the proof of the theorem, we will apply the dominated 
convergence theorem. We have the following estimate for y < St- 



OR 

M Y e , a <Pe(y ~ a) — (Y £ , a ,T y )4){a) da 



< 



oo SUp 

s,t<T 



OR, , 

m 1 ^ 



sup|Mj|, 

t<T 



which allows us to commute the limit (3.10) with the integral with respect 
to the measure P x dr y on the set {(a;, y) :y < St(u>)}. In this way we get 

E{M Ty )^{y)dy 
ip(y) dy 

( [ St f 1 OR \ 

~ XE [j ^(f) J M ZT y+ +(l-z)Ty-fc( ZT y+ + (1 - Z)Ty,T y )dzdT y \. 

Approximating l[o )0 ] by a sequence of smooth functions (tp n ,n > 1) and 
letting T tend to infinity completes the proof. □ 
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If we assume that the partial derivative Q^(t,s) is nonnegative, then we 
can derive the following result. 

Proposition 3.5. Assume that X satisfies hypotheses (HI), (H2) and 
(H3). If Qfc{s,t) > 0, then, for all a, a > 0, we have 

(3.11) E{exp{-aV Ta ))<e~ aV **. 

PROOF. Since ^(t,s) > 0, we obtain 

E(M Ta )<l, 

that is, 

£(exp(Aa-iA 2 Kj)<l, 

or 

E(exp(-aV Ta )) < e~ aV ^. 

The result follows. □ 

The above proposition means that the Laplace transform of the random 
variable V Ta is dominated by the Laplace transform of r a , where r a is the 
hitting time of the level a for the ordinary Brownian motion. This domi- 
nation implies some consequences on the moments of V Ta . In fact, for any 
r > 0, we have, multiplying (3.11) by a r , 

1 r°° 

E{v ^ r) = W)l E{e " Ta )<* r ~ Xda 

1 roo 

(3.12) e-^a-Ua 

_2T(r+l/2) (j _ 2r 



/7T 

On the other hand, for < r < 1, 



E(V T r J = — — -J (l-B(e-^))a-Ma 

(3.13) ° 

r r°° i — 

In particular, for r £ [1/2, 1), = +°°- 



Remark 3.6. If X is the standard Brownian motion, its covariance 
s At does not satisfy condition (HI), but we still can apply our approach. 
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It is known from [3] that dr a has no absolutely continuous part and that 
{a,T a = t+} is a Cantor set, hence, of zero Lebesgue measure. It follows 
from this observation and from (3.10) that 

J E{M Ty )^{y)dy = J ^(y)dy. 

Choosing -0 = l[o,a] yields to the expected result: 

E ([ eXy ~ {x2 ' 2)v{Ty)d y 

If X has independent increments and satisfies (H3), then 

£ ( e -(AV2)V(r ) ) = e -A« 

This follows easily from the fact that X can be written as a time-changed 
Brownian motion. 




Remark 3.7. Consider that X is a fractional Brownian motion of Hurst 
index H = 1. Then R(s, t) = st, and consequently, Xt = Yt, where Y is a one- 
dimensional standard Gaussian random variable. Then, r a = r a + = a/Y + . It 
is then easy to compute the Laplace transform of r a and we obtain 

(3.14) £(exp(-ar a 2 )) = \e~ a ^. 

We show now that our formula also yields to the right answer. We just note 
that {y i — > T y ) is continuous. This entails that 

9R . , . . dR . . l T w/ x 

— (zT y+ + (l-z)T y ,Ty) = -g^( T y^y) = 2 V W 



and 



(3.15) 



Set 



then 



= a - ±e(£ exp (\y - ^ K(r„)) V'(r y ) dr y 



(a, \) = E[ exp ( Aa - — F(v 



,3,6) » M) 



12 



L. DECREUSEFOND AND D. NUALART 



Substitute (3.15) into (3.16) to obtain 

-— = 2A<I> - A. 

da 

Then, there exists a function C(A) such that 

¥(a, A) = - + C(A)e 2Aa so that ^fexp (-— tA) = -e~ Xa + C(A)e Aa . 



o 



2 v ' V V 2 a J J 2 

By dominated convergence, it is clear that, for any A, 

thus, C(A) = and 

B K-T*))=k*- 

Changing A 2 /2 into a gives (3.14). 

Remark 3.8. Consider the case of a fractional Brownian motion with 
Hurst parameter H > 5. Conditions (HI), (H2) and (H3) are satisfied and 
we obtain 

° a E{M Ty )dy 

= a- \HE(J a J o M ZTy++{1 _ z)Ty ([zT y+ + (1 - z)r y ] 2H - 1 

-\z{j y + ~ Ty)\ 2H ~ 1 )dzdTy 

Moreover, E(e~ aT " H ) < e _av/2 ^, and therefore, E(t%) < 00 if p < H. Accord- 
ing to (3.13), E(t%) is infinite if pH > 1/4 and (3.12) entails that r a has 
finite negative moments of all orders. 
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